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Abstract 

Background: Environmental justice research has shown that nnany connnnunities of color and low-income persons 
are differentially burdened by noxious land uses including Toxic Release Inventory (TRI) facilities. However, limited 
work has been performed to assess how these populations tend to be both overburdened and medically 
underserved. We explored this "double disparity" for the first time in Maryland. 

Methods: We assessed spatial disparities in the distribution of TRI facilities in Maryland across varying levels of 
sociodemographic composition using 2010 US Census Health Professional Shortage Area (HPSA) data. Univariate 
and multivariate regression in addition to geographic information systems (CIS) were used to examine relationships 
between sociodemographic measures and location of TRI facilities. Buffer analysis was also used to assess spatial 
disparities. Four buffer categories included: 1) census tracts hosting one or more TRI facilities; 2) tracts located more 
than 0 and up to 0.5 km from the closest TRI facility; 3) tracts located more than 0.5 km and up to 1 km from a TRI 
facility; and 4) tracts located more than 1 km and up to 5 km from a TRI facility. 

Results: We found that tracts with higher proportions of non-white residents and people living in poverty were 
more likely to be closer to TRI facilities. A significant increase in income was observed with an increase in distance 
between a census tract and the closest TRI facility. In general, percent non-white was higher in HPSA tracts that 
host at least one TRI facility than in non-HPSA tracts that host at least one TRI facility. Additionally, percent poverty, 
unemployment, less than high school education, and homes built pre-1950 were higher in HPSA tracts hosting TRI 
facilities than in non-HPSA tracts hosting TRI facilities. 

Conclusions: We found that people of color and low-income groups are differentially burdened by TRI facilities in 
Maryland. We also found that both low-income groups and persons without a high school education are both 
overburdened and medically underserved. The results of this study provide insight into how state agencies can 
better address the double disparity of disproportionate environmental hazards and limited access to health care 
resources facing vulnerable communities in Maryland. 
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Introduction 

Environmental injustice is driven by privilege, power- par- 
ticularly structural and environmental racism which are 
embedded in our regulatory schema, zoning, planning and 
community development processes [1,2]. Additionally, un- 
healthy geographies that concentrate environmental, so- 
cial, and health risks in urban and rural areas are 
produced and are known as riskscapes' [3,4]. The original 
Toxic Waste and Race in America report published in 
1987, was the first report to demonstrate that many eco- 
nomically underserved populations and people of color 
communities are disproportionately impacted by locally 
unwanted land uses (LULUs) [5]. Since the 1987 report 
and the recent twenty-year anniversary report [5], re- 
searchers have shown that these disparities persist, with 
low-income persons and populations of color continuing 
to live in communities with a differential burden of LULUs 
including toxic release inventory (TRI) facilities [6-18], 
landfills [5], incinerators [5], hazardous waste sites [5,18], 
sewer and water infrastructure including sewer and water 
treatment plants [7,8,19,20], coal-fired plants [5], industrial 
animal operations [21,22], and Superfund sites [23-25]. 
This disproportionate burden can lead to increased expos- 
ure to harmful environmental conditions and chemical, 
physical, and biological agents for impacted communities 
[1,2,26-28]. 

Previous research has also shown that populations of 
color and low-income groups living in poor environmental 
conditions have health risks due in part to various social 
determinants of health including segregation, racism, so- 
cioeconomic status (SES), income inequality, and inequities 
in planning and zoning [1-4,6,14,19,20,28-31]. Studies have 
shown that underlying social and economic vulnerabilities 
contribute to increased health disparities [29,31,32], which 
further enhance the long-term effects of environmental 
injustice. Environmental justice communities are also af- 
fected by a higher concentration of psychosocial stress 
[1,28,29] that can lead to an increase in community-level 
and individual-level stress. 

A potential environmental justice issue in the state of 
Maryland is the distribution and concentration of TRI fa- 
cilities. Previous studies in New York [7], South Carolina 
[9,10], Oregon [11,12], California [13,18], and the entire 
United States [6,8,15] have demonstrated the dispropor- 
tionate burden of TRI facilities in low-income and non- 
white communities (often using census tracts as the unit 
of analysis). Ringquist found that TRI facilities were found 
in zip codes with large populations of people of color [15]. 
Racial composition of neighborhood was found to be a 
stronger contributor to the trend associated with the dis- 
tribution of environmental risk than class [15]. Abel found 
that people of color and low-income residents were 
disproportionately closer to TRI facilities in metropolitan 
St. Louis [16]. Spatial concentration of residents of color 



averaged nearly 40% within one km of St. Louis TRI sites 
compared to 25% in other locations [16]. Pricker and 
Hengarter report that the racial/ethnic composition of a 
census tract in Metropolitan New York was positively as- 
sociated with the presence of LULUs including TRI sites 
[7]. In metropolitan New York, the Hispanic population 
resided in neighborhoods closer to undesirable sites than 
other racial/ethnic groups [7]. Using 1990 US Census and 
1990 TRI data, Daniels and Friedman observed a positive 
relationship between proportion Black residents and toxic 
releases to air [17]. 

Wilson et al. found significant burden disparities [9], 
where more TRI facilities were located in census tracts 
with higher non-white and low-income populations for the 
state of South Carolina and Metropolitan Charleston. In 
addition to this work, other researchers have documented 
similar racial and income disparities among communi- 
ties hosting TRI facilities [11,12]. Neumann et al discov- 
ered that TRI facilities were located disproportionately 
in people of color neighborhoods and in areas with 
lower incomes compared to those in the surrounding 
counties [11]. 

Miranda et al. took these analyses one step further to 
scrutinize the effects of new TRI reporting requirements 
implemented in December 2006 which reduced reporting 
requirements for certain chemicals released in limited 
quantities [33]. Specifically, the study found that facilities 
given permissions to use a short reporting form were dis- 
proportionately located in majority non- white census 
tracts [33]. As a result, residents were losing access to sali- 
ent information regarding chemical releases compared to 
their white counterparts who were mostly located in areas 
with more stringent reporting requirements. 

The paucity of information on chemical releases is prob- 
lematic because many of the chemicals typically emitted 
by TRI facilities and other LULUs have been linked to ad- 
verse cancer and non-cancer health effects [34,35] which 
may significantly impact people who live near these envir- 
onmental hazards. For example, emissions from TRI facil- 
ities may include harmful substances such as benzene, 
cadmium, toluene, and mercury, among other chemicals 
[36,37]. Moreover, populations exposed to TRI-related 
chemicals may have an increased risk of adverse health 
outcomes such as low birth weight, asthma, and cancer 
[13,38-40] estimated by previous research using risk as- 
sessment methodology. 

In addition to being differentially burdened by environ- 
mental hazards and LULUs, the lack of access to saluto- 
genic infrastructure (e.g., positive and health-promoting 
features of the built and social environment) [1,2,41] is a 
major environmental justice issue for people of color 
communities and economically disadvantaged populations. 
Many persons of color live in socially disadvantaged areas 
with limited access to primary care resources [42-46]. 
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Access to adequate healthcare may be a major problem for 
economically disadvantaged in communities with pre- 
existing burden, exposure, and environmental health dis- 
parities. A growing body of literature has examined the 
distribution of health enriching resources and medical care 
services across varying racial/ethnic and SES composition 
at the neighborhood level Limited access to hospitals and 
medical professionals and lower quality of care both play 
major roles in health outcomes and disparities in disadvan- 
taged neighborhoods [47-52]. Being both disadvantaged 
and medically underserved means disadvantaged popula- 
tions may have higher rates of chronic conditions, more 
drug use, emotional problems, and worse health behaviors 
than other populations [53], but lack of access to compe- 
tent high quality care may mean lower immunizations 
rates for children [54] and more hospitalizations for treat- 
able and preventable conditions [55]. Taken as a whole, 
the differential burden of pollutogens and access to saluto- 
genic resources has important implications for health and 
environmental health disparities [1,2,41]. 

The issue of how people of color neighborhoods and 
low-income populations are disproportionately burdened 
by LULUs such as TRI facilities and are potentially under- 
served due to poor access to medical infrastructure is im- 
portant for a state like Maryland with a number of racial/ 
ethnic and SES -related health disparities. The 2010 MD 
Plan to Eliminate Minority Health Disparities (MPEMHD) 
lists racial/ethnic disparities in healthcare utilization, ac- 
cess to primary care, and the burden of all-cause mortality, 
heart disease, renal disease, hypertension, obesity, HIV/ 
AIDS, and asthma as critical areas for improvement [56]. 
Three of the four counties in MD with the highest popu- 
lation of persons of color (Baltimore City, Charles County, 
Montgomery County, and Prince Georges County) all 
have more than 50% non-white residents and the high- 
est total environmental releases (Charles County, Prince 
Georges, and Baltimore City) [56]. As further proof that 
these health disparities are a serious problem, the Maryland 
Health Improvement and Health Disparities Reduction Act 
(MHIHDRA) was recently passed to address the aforemen- 
tioned disparities [56]. 

In addition, the MD Department of Health and Mental 
Hygiene (MDHMH) uses its Environmental Public Health 
Tracking Network (EPHTN) to offer web-based data shar- 
ing tools for residents, policymakers, and other public offi- 
cials to create their own maps and charts documenting 
environmental health disparities in their service areas 
[57,58]. Information available for analysis through the MD 
EPHTN includes the following: 1) childhood blood level 
testing, 2) myocardial infarction and asthma-related hos- 
pitalization data, and 3) low birth weight from state birth 
certificate records [57] . Thus, running a query reveals sig- 
nificant disparities in blood lead levels among counties; 
with Baltimore City having a much higher rate than other 



counties (449 1-year olds in 2008 had elevated blood lead 
levels versus many other counties that had none) [59]. An- 
other query revealed that asthma hospitalization discharge 
rates were higher among blacks compared to whites across 
the state (36.79 per 10,000 vs. 11.20 per 10,000, respect- 
ively), and this disparity was more pronounced in certain 
areas (32.33 per 10,000 vs. 2.88 per 10,000 in Baltimore 
City and 10.31 per 10,000 vs. 1.00 per 10,000) [56,59]. 

The purpose of this study was to assess whether TRI 
facilities in Maryland were more likely to be located in 
census tracts with higher proportions of black, non- 
white, low-income, or less educated persons. In addition, 
we assessed whether populations near TRI facilities had 
limited access to health care infrastructure as indicated 
by health professional shortage area (HPSA) designation 
at the census tract level. By assessing both presence of 
TRI facilities and HPSA designation, we assessed the po- 
tential "double disparity" of being environmentally over- 
burdened and underserved in terms of health care 
access across neighborhoods with varying sociodemo- 
graphic composition. 

Methods and materials 

Study area 

The state of Maryland (MD) is ranked 42^*^ in size among 
states in the U.S., but 19^^ in population, which makes it 
one of the more densely populated states in the nation 
[60]. The population within MD is concentrated in two 
main areas: 1) around the harbor in Baltimore County and 
Baltimore City and 2) Montgomery County and Prince 
Georges County, MD near Washington, DC. According to 
the 2010 US Census, there were 5,773,552 people living in 
MD with 61.1% white and 30.0% black [60]. Furthermore, 
populations of color living in MD are highly concentrated 
in these two areas. Baltimore City is 72% non-white and 
Prince Georges County is 85% non-white, while MD is 
39% non-white as a whole [60]. 

The number of people living in poverty is also un- 
evenly distributed throughout the state. For example, 
Maryland ranks 3"^*^ in the nation in median household 
(HH) income with only 8.6% of the state living in pov- 
erty (compared to 13.8% nationwide) [60]. In Baltimore 
City, the most densely populated urban area in the state, 
21.3% of residents live below the federal poverty line 
[60]. In two other distinct parts of the state, Allegany 
County in western MD has 14.5% of its residents living 
in poverty and Dorchester County on the eastern shore 
has 13.4%, both well above the state poverty rate [60]. 

Sociodemographic (SOD) Measures 

This study used key demographics for MD modeled as 
quartiles from 2010 census data [9,61,62]. While SOD 
information is available at various geographic scales (ZIP 
code tabulation areas (ZCTAs), tracts, block groups, and 
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blocks), we utilized census data at the tract level to enu- 
merate the following population characteristics: race/ 
ethnicity (% non-white includes all other races including 
Hispanics except non-Hispanic white; and % Hispanic) 
and variables related to socioeconomic status (SES). SOD 
measures included in our study were poverty (% popula- 
tion below poverty line), education (% population age 
greater than 25 years with < high school (HS) education), 
unemployment (% of population 16 years and older who 
were unemployed), homeownership (% of homes occupied 
by owners), and homes built before 1950 (% of homes 
built pre- 1950). The variable related to income was me- 
dian HH income. Median HH income, % poverty, house 
construction year, % unemployment and educational at- 
tainment (i.e., % without a HS diploma) were calculated 
using the 2006 to 2010 American Community Survey 
(ACS) 5 -year estimates. 

USA today diversity index 

The USA Today Diversity equation measures the distribu- 
tion of multiple races. Specifically, it calculates the prob- 
ability that any two people randomly selected in an area 
are from different races or ethnic groups [63]: 

USA Today Diversity = l-((l-qh)^ + x ^p^ 

where qh is the % of Hispanics in an area and Pi includes 
% White, African-American, Asian, American Indian 
and Alaska Native, and Native Hawaiian and other Pa- 
cific Islander. The higher the Diversity Index, the more 
diverse a population in a particular region ranged from 0 
to 100. The Diversity Index is well-suited for character- 
izing racial and ethnic diversity at the tract level, while 
other indices are more appropriate for examining segre- 
gation at larger geographic scales, such as the MSA. In 
this study, the Diversity Index was calculated at the cen- 
sus tract level. 

Toxics release inventory (TRI) 

The TRI database was established by Section 313 of the 
1986 Emergency Planning and Community Right- to- 
Know Act (EPCRA) [64]. Estimations of the mass of 
disposal or other release of over 650 chemicals were re- 
ported to the TRI by each facility (2010 Toxics Release 
Inventory national analysis overview). The 2010 single 
Facility Registry System (FRS) state files were down- 
loaded from the USEPA and TRI facilities were selected 
from the FRS which contains their respective latitude 
and longitude coordinates. Distance between a TRI facil- 
ity and nearest census tract was calculated in ArcGIS 10 
(esri, Redlands, CA). A TRI facility located in a census 
tract or on the boundary of the census tract was 
assigned a distance of 0 which means that the census 
tract 'hosts' the TRI facility. 



Healthcare infrastructure 

The Department of Health and Human Services (DHHS) 
created Health Professional Shortage Area (HPSA) desig- 
nation to identify areas facing a critical shortage of pro- 
viders [42-46,65]. A HPSA can be a distinct geographic 
area (such as a county), a specific population group within 
an area (such as low-income individuals), or a specific 
health care facility [66]. We obtained 2010 HPSA data for 
the state of Maryland and categorized each census tract in 
Maryland either as a HPSA census tract or non-HPSA 
census tract. 

Statistical and geographic methods 

To assess the proximity from TRI facilities to different 
sociodemographic features, we used two approaches. In 
the first approach, census tracts were grouped by their 
distance to the nearest TRI facility. A distance was mea- 
sured along a straight-line path between a facility and 
the point closest to the facility on the boundary of a cen- 
sus tract. Based on this definition, we grouped census 
tracts into four groups (or four distance bands), defined 
as: band 1 composed of census tracts hosting one or 
more TRI facilities (at least one TRI in the census tract 
and the distance was 0), band 2 included census tracts 
whose distance to the closest TRI facility was greater 
than 0 and up to 0.5 km, band 3 greater than 0.5 km 
and up to 1 km, and finally band 4 greater than 1 km 
and up to 5 km. Census tracts whose distance to the 
nearest TRI facility was greater than 5 km but within 
10 km were excluded from this analysis due to the small 
number of census tracts in this distance band. Then 
within each band, the mean percentage of each SOD 
measure for those census tracts was calculated. This task 
was performed to evaluate how SOD composition chan- 
ged with the change in distance. When disproportional- 
ity did not exist with respect to a specific SOD group, 
we expected that the mean SOD measure (by percent- 
age) in each distance band would correspond to the 
mean percentage calculated for the entire state. Student 
t-tests were used to test this hypothesis of equality. 

In the second approach, census tracts were grouped 
based on the percent number for a specific SOD vari- 
able. This was done to evaluate how distance changed 
with respect to the change in the composition of popula- 
tion for a specific SOD group. To this end, we divided 
the census tracts into four groups (Ql to Q4) using 
quartiles of a SOD variable measured by percentage 
across the state of Maryland. Then in each group of cen- 
sus tracts, the average distance between each tract and 
the closest TRI facility was calculated. 

To quantify the relationship between the distribution 
of TRI facilities and SOD factors at the census tract 
level, we first applied univariate linear regression to test 
the association of distance from census tracts to the 
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closest TRI facility (dependent variable) with each indi- 
vidual SOD factor (independent variable). This was 
then followed by a multivariable linear regression 
model with all SOD factors included in order to each 
SOD factors effect after adjusting for other SOD fac- 
tors, i.e., yi = |3o+ ^^X + 8i where yi denotes distance, X 
denotes a vector of SOD factors, the ^ vector is for the 
SOD factor effects, and e is the random error. To elim- 
inate redundant variables, a stepwise variable selection 
approach based on the Bayesian Information Criterion 
(BIC) was applied to finally select the most important 
variables. 

We hypothesized that low-income and people of color 
communities are both overburdened by environmental 
hazards and unhealthy land uses and underserved by 
health infrastructure known as the "double disparity". 
To test the difference in SOD composition between 
HPSA tracts and non-HPSA tracts, we compared the 
average SOD measures in HPSA tracts and non-HPSA 
tracts in two distance band areas, one area covered cen- 
sus tracts hosting a TRI facility and the other area in- 
cluding census tracts within 1 km to 5 km to the nearest 
TRI facility. For each area, a student t-test was used to 
test the difference in composition for each SOD measure 
between HPSA tracts and non-HPSA tracts. 

All of the above calculations and hypothesis testing were 
performed using R version 2.15.0 [11,12,67]. Statistical 



significance level was set at 0.05. TRI facilities were 
mapped and overlaid by select SOD features (% non- 
white and % poverty) in ArcGIS 10 (esri, Redlands, CA). 
Choropleth maps were created to illustrate the spatial 
relationship between TRI facilities and sociodemo- 
graphic composition using quartiles. 

Results 

Figures 1 and 2 were created to show the spatial distri- 
bution of TRI facilities in relation to composition of 
various sociodemographic groups at the census tract 
level. There were 525 TRI facilities located in 259 census 
tracts in MD out of a total of 1390 census tracts. In 
addition, there was one census tract located close to the 
Port in South Baltimore that hosts 14 TRI facilities 
which was the largest number of TRI facilities found in 
one tract in the state. The choropleth maps show clus- 
ters of TRI facilities in the Baltimore Metropolitan Stat- 
istical Area (MSA), Washington County (western MD), 
and Wicomico County (Eastern Shore). 

For the purpose of this analysis, we only present map- 
ping results for % non-white and % poverty due to the 
fact that previous studies have shown a positive relation- 
ship between number of TRI facilities and presence of 
non-whites and low-income persons. For % non-white, 
137 TRI facilities were located in the first quartile 
followed by 165, 164, and 47 facilities in the second. 
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Figure 1 Choropleth Map of TRI Facilities in Maryland by Quartiles for Percent Non-White (2010 US Census). 
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Figure 2 Choropleth map of TRI facilities in Maryland by quartiles for percent poverty (2010 US census). 



third, and fourth quartiles, respectively. There were 12 
TRI facilities not included on the maps due to a loca- 
tional error or they were located in census tracts where 
no people lived. The census tracts in the fourth quartile 
had the largest non-white population but the fewest TRI 
facilities compared to other quartiles. Regarding % living 
in poverty, 70 TRI facilities were located in the first 
quartile followed by 108, 157, and 178 facilities in the 
second, third, and fourth quartiles, respectively. Both 
figures show evidence of differences in the concentration 



of TRI facilities across varying levels of sociodemo- 
graphic composition at the census tract level. Figure 2 
indicated a clear linear relationship between % living in 
poverty and TRI facilities. Thus, as % persons in poverty 
increased, the number of TRI facilities increased. 

Statistical summaries and tests indicated that the popu- 
lation distribution for different SOD variables changed 
with an increase in distance to TRI facilities (Table 1). For 
example, mean % non-white in census tracts having at 
least one TRI facility in their geographic boundary (host) 



Table 1 Mean distribution of sociodemographic measures by TRI facility buffer zones in Maryland (2010 census) 



Sociodemographic measure 


State 


Host (a census tract that hosts at least one TRI facility) 


0-0.5 km 


0.5-1.0 km 


1-5.0 km 


% Hispanic 


7.8 


6.9 


8.5 


7.8 


8.6 


% Non-white 


46.1 


38.4** 


49.2 


54.8** 


49.1 


% Poverty 


9.6 


ir 


10.8 


12.7** 


8.3* 


% Unemployment 


7.1 


7.2 


8** 


8.3** 


6.6* 


% < HS education 


13.2 


15.3** 


14.6 


16** 


1 1 .7** 


% Homeownership 


67.1 


65.3 


63.3* 


59.3** 


68.9 


% Homes built pre-1950 


20.8 


25.4** 


24.6* 


28.7** 


16.9** 


^Diversity index 


0.43 


0.43 


0.46* 


0.43 


0.44 


Median HH income 


74810 


65239** 


69399* 


64362** 


81190** 


MD census tracts (N) 


1390 


259 


218 


192 


613 



*p-value <0.05 ^^value <0.01. 

^The Diversity index was 0.428 for state and 0.433 for 0.5 to 1.0 km group. 
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(38.4%) was significantly lower than the mean % non- 
white in all census tracts in Maryland. With an increase in 
distance to the nearest TRI facility, % non-white increased 
and becomes statistically significantly higher than the 
statewide mean for % non-white. The census tracts with 
distance to the nearest TRI facility in the third distance 
band (>0.5 km and <1 km) had the highest % non-white 
(54.8%) which was 8.7% higher than the statewide average 
Results indicate that on average, the highest proportion of 
non-white residents was in tracts that were between 0.5- 
1 km distance from a TRI facility. However, there was no 
statistically significant difference between statewide mean 
% Hispanic and mean % Hispanic in host tracts or other 
tracts grouped by distance. 

Percent living in poverty in host tracts and % poverty 
of census tracts in the second distance band was almost 
the same (11% and 10.8%, respectively). They both were 
higher than the statewide average (9.6%). Percent living 
in poverty increased to 12.7% in the third distance band 
(>0.5 km and <1 km) and decreased to 8.3% for the 
fourth distance band (>1 km and <5 km). In addition, 
these results reveal that there were more persons in pov- 
erty living in areas closer to TRI facilities especially in 
census tracts within the third distance band (>0.5 km 
and <1 km). We observed a similar pattern for changes 
in % less than HS education as distance from census 
tracts to the nearest TRI facility increased (Table 1). For 
% unemployment, there were no statistically significant 
differences between statewide average percentage and in 
areas hosting TRI facilities (7.1% for statewide average 
and 7.2% for host tracts). 

For housing-related variables, we observed lower % 
homeownership in host tracts compared to the statewide 
average, but the difference was not statistically signifi- 
cant. As the distance from census tracts to the nearest 
TRI facility increased, mean % homeownership de- 
creased from 63.3% to 59.3% from the second distance 
band (>0 km and <0.5 km) to the third distance band 
(>0.5 km and <1 km) both of which were statistically 
lower than the statewide average (67.1%). After 1 km, % 



homeownership increased to 68.9% which was 1.8% 
higher than the statewide average. The changes in % 
homes built pre- 1950 were similar to changes in % pov- 
erty (Table 1). 

As for the Diversity index, no clear pattern was ob- 
served. Statistically significant differences were observed 
in the 0-0.5 km buffer (0.46). At the census tract level, 
the statewide average median HH income ($74,810) was 
almost $10,000 higher than median HH income for cen- 
sus tracts hosting a TRI facility and in areas in the sec- 
ond distance band (>0.5 km and <1 km). The median 
HH income increased to $81,190 in areas in the fourth 
distance band (>1 km and <5 km). These results indicate 
that TRI facilities possibly cluster in low-income areas 
that host the LULU or in areas located at least 1 km 
away from the nearest TRI facility. 

Table 2 shows the mean distance between TRI facil- 
ities and each group of census tracts defined by quartiles 
of SOD measures. When considering race, the mean dis- 
tance from TRI facilities to census tracts grouped by % 
Hispanic from Ql to Q4 decreased 20%, while % non- 
white from Ql to Q4 decreased 40% (2.5 km vs 1.5 km). 
These results indicate that census tracts with a higher % 
non-white population were located closer to TRI facil- 
ities. For % poverty, the mean distance from census 
tracts to the nearest TRI facility decreased 50% from Ql 
to Q4. A similar pattern in distance change across differ- 
ent quartile groups of census tracts was observed for % 
unemployment and % less than HS education (Table 2). 
As for % homeownership, the higher the percentage, the 
farther the distance to the nearest TRI facility (mean 
distance of 1.2 km for Ql and 2.7 km for Q4). Percent 
homes built pre-1950 showed a non-linear pattern. The 
distance first increased from Ql to Q2, and then de- 
creased from Q3 to Q4. The Diversity index showed a 
decreasing pattern of distance from Ql to Q4, consistent 
with the pattern revealed by % poverty, % unemploy- 
ment, and % less than HS education. The mean distance 
to the nearest TRI facility for census tracts grouped by 
quartiles of median HH income showed a linear increase 



Table 2 Mean distance to TRI facilities by quartiles for various sociodemographic measures in Maryland (2010 census) 



Sociodemographic measures 


Ql (km) 


Q2 (km) 


Q3 (km) 


Q4 (km) 


Q4/Q1 


% Hispanic 


2.1 


1.9 


1.6 


1.6 


0.8 


% Non-white 


2.5 


1.7 


1.4 


1.5 


0.6 


% Poverty 


2.3 


1.9 


1.8 


1.1 


0.5 


% Unemployment 


2.2 


1.9 


1.8 


1.3 


0.6 


% < HS education 


2.3 


2.2 


1.7 


1 


0.4 


% Homeownership 


1.2 


1.3 


2 


2.7 


2.3 


% Homes built pre-1950 


2 


2.2 


1.8 


1.1 


0.6 


Diversity index 


2.2 


2 


1.4 


1.4 


0.6 


Median HH income 


1.1 


1.5 


1.9 


2.6 


2.4 



Wilson et at. Environmental Health 2014, 13:26 
http://www.ehjournal.net/content/1 3/1/26 



Page 8 of 1 2 



Table 3 Linear regression of decay in distance to TRI facilities by exposure factors in Maryland (2010) 





Univariate linear regression 


Multivariate linear regression 


Sociodemographic measures 


Beta coefficient 


p-value 


Beta coefficient 


p-value 


% Hispanic 


-0.014 


0.01 


0.028 


<0.001 


% Non-white 


-0.01 1 


<0.001 






% Poverty 


-0.043 


<0.001 






% Unemployment 


-0.071 


<0.001 






% < HS education 


-0.044 


<0.001 






% Homeownership 


0.022 


<0.001 


0.006 


<0.046 


% Homes built pre-1950 


-0.019 


<0.001 


-0.017 


<0.001 


Diversity index 


-1.409 


<0.001 


-2.523 


<0.001 


Median HH income 


1.15X 10"^ 


<0.001 


0.85 X 10"^ 


<0.001 


HPSA 


-0.717 


<0.001 







-: variables in the Initial multivariate model before stepwise model selection (both direction) based on BIC and confounder selection. 



from Ql to Q3 with a rate of 0.4 km. The mean distance 
then increased to 2.6 km in the Q4 group which was al- 
most 2.4 times farther away than the Ql census tract 
group. 

In the univariate regression model, all SOD mea- 
sures were statistically significantly associated with 
distance to TRI facilities (Table 3). Across the state of 
MD, census tracts with a higher percentage of non-white 
residents demonstrated resistance to decaying distance 
between adjacent TRI facilities (Beta coefficient = -0.011; 
p < 0.001). We observed similar effects in the same direc- 
tion for % poverty, % unemployment, % < HS education, % 
homes built pre-1950 and Diversity Index (Beta coeffi- 
cient = -0.043, -0.071, -0.044, -0.019, -1.409, respectively, 
with p < 0.001 in all these tests). Conversely, 1% increase 
in homeownership in census tracts would increase the dis- 
tance by 0.022 km (p < 0.001). For income-related variables, 
a significant increase in income was observed with an in- 
crease in distance between a census tract and the closest 
TRI facility (1.15 xlO'^ for median HH income, p < 0.001). 

In the multivariate regression model, after performing 
stepwise variable selection (Additional file 1: Table SI) 
and evaluating the confounding effects of the SOD fac- 
tors (Additional file 2: Table S2), % non-white, % homes 
built pre-1950. Diversity index and median HH income 
were shown to be significantly associated with TRI facil- 
ity distance from census tracts (p < 0.001) with % home- 
ownership as a potential confounder. The direction of 
the effects remained the same except for % Hispanic. 
The direction of the association between distance to TRI 
facilities and % Hispanic changed from negative to posi- 
tive. This is likely due to the adjustment for potential 
confounders (e.g., % unemployment, % less than HS edu- 
cation, and % homeownership). 

Table 4 presents mean distribution of SOD measures 
in HPSA and non-HPSA tracts based on distance to TRI 



facilities from census tracts. There were 57 HPSA tracts 
and 202 non-HPSA tracts that host TRI facilities. There 
were 117 HPSA tracts and 496 non-HPSA tracts in areas 
where the distance to the nearest TRI facility from a 
census tract was from 1 km to 5 km. In general, % non- 
white was higher in HPSA tracts that hosted at least one 
TRI facility than in non-HPSA tracts that hosted at least 
one TRI facility. Additionally, % poverty, % unemploy- 
ment, % less than HS education, % homeownership, and 
% homes built pre-1950 were higher in HPSA tracts 
hosting TRI facilities than in non-HPSA tracts hosting 
TRI facilities. All the means of SOD measures in HPSA 
and non-HPSA tracts were statistically different except 
for % Hispanic in both host tracts and tracts with near- 
est TRI facilities at a distance from 1 km to 5 km away. 

Overall, statistically significant higher levels of % non- 
white, % poverty, % unemployment, % less than HS 



Table 4 Mean distribution of sociodemographic measures 
by TRI facility buffer zones for 2010 Maryland HSPA 
tracts and non-HPSA tracts 



Sociodemographic measures 


Host 


1 km - 5 km 
buffer 




HPSA/Non-HPSA 


HPSA/Non-HPSA 


# Census tracts 


57/202 


1 1 7/496 


% Hispanic 


8.3/6.6* 


10.3/8.2* 


% Non-white 


55.8/33.4 


79.8/41.9 


% Poverty 


18.8/8.8 


16/6.5 


% Unemployment 


11.2/6.2 


11.1/5.5 


% Less than HS education 


22.2/13.3 


21/9.5 


% Homeownership 


53.4/68.6 


47.9/73.8 


% Homes built pre-1950 


41 .5/20.9 


27.6/14.4 


Diversity index 


0.41/0.43* 


0.37/0.45 


Median HH income 


47428/70202 


48723/88798 



^Statistically insignificant at the level of 0.05. 
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education and % homes built pre- 1950 were observed 
in HPSA tracts than in non-HPSA tracts, regardless 
of whether or not those areas hosted a TRI facility. 
Additionally, % homeownership was lower in HPSA tracts 
compared to non-HPSA tracts. For the Diversity Index, in 
areas within 1 km to 5 km to a TRI facility, non-HPSA 
tracts had a higher index than HPSA tracts. Median HH 
income was higher in HPSA tracts than non-HPSA tracts, 
regardless of whether or not the tract hosted a TRI facility. 

Discussion 

Our results primarily indicate that people of color, low- 
income populations, and persons with less than HS educa- 
tion are located closer to TRI facilities than other groups 
or there are greater numbers of people of color and low- 
income persons in census tracts in areas that host TRI 
facilities. However, we did observe a lag effect of TRI facil- 
ities on the distribution of non-whites across different dis- 
tance bands with higher percent non-white in the third 
distance band (>0.5 km and <1 km) compared to the host 
band (0 km) and the second distance band (>0 km 
and <0.5 km). For % poverty and % less than HS education, 
we did not observe this lag effect. Univariate regression 
results reveal statistically significant inverse relationships 
between distance to TRI facilities and % Hispanic, % non- 
white, % poverty, % unemployment, % less than HS edu- 
cation, % homes built pre-1950, and Diversity index (all 
statistically significant). Conversely, a positive relationship 
was observed between distance to TRI facilities and % 
homeownership and median HH income. In the multivari- 
ate models, results reveal statistically significant inverse re- 
lationships between distance to TRI facilities and % homes 
built pre-1950 and Diversity Index, and a positive relation- 
ship between distance and median HH income and % 
homeownership. Percent Hispanic changed from negative 
to positive which may be due to adjustment of other 
SOD factors in the multivariate model that could act as 
confounders. 

The results of our study mirror results of previous re- 
search that demonstrated burden disparities in the dis- 
tribution of TRI facilities and other LULUs based on 
race/ethnicity and class [6-18]. The results of this study 
are most similar to results of recent work in metropol- 
itan Charleston [9] and St. Louis [16] where more TRI 
facilities were located in census tracts with higher non- 
white and low-income populations. Additionally, Fricker 
and Hengarter found a direct relationship between racial 
composition and presence of TRI facilities [7]. Both 
Ringquist and Neumann et al found TRI facilities were 
located in people of color neighborhoods [11,15]. Unlike 
Ringquist [15], we did not use zip codes but census 
tracts as the unit of analysis. 

A benefit of including HPSA data is that we can com- 
pare HPSA census tracts to non-HPSA census tracts to 



assess relative differences in the mean distribution of 
SOD measures by TRI host and buffer zones. This ap- 
proach is to understand if people of color and disadvan- 
taged groups are both overburdened and medically 
underserved which is a problem when addressing envir- 
onmental health disparities. For example, we observed 
that % poverty in HPSA census tracts (tracts with insuf- 
ficient access to primary care) hosting a TRI facility was 
more than twice the % poverty in the equivalent non- 
HPSA census tracts (18.8% and 8.8%, respectively). Fur- 
thermore, we found that % < HS education population in 
HPSA census tracts hosting TRI facilities was almost 
twice that of non-HPSA census tracts (22.2% and 13.3%, 
respectively). 

Aside from disparities in mean poverty between HPSA 
and non-HPSA TRI host census tracts, non-white popu- 
lations had a higher percentage of persons hosting TRI 
facilities in HPSA versus non-HPSA tracts (non-white: 
55.8%, 33.4%, respectively). Additionally, we observed 
disparities in % non-white, % poverty, % unemployment, 
and % less than HS education and also median HH in- 
come in the 1 km to 5 km buffer for HPSA vs non- 
HPSA tracts. The disparity in race, SES, and educational 
attainment in HPSA versus non-HPSA census tracts is 
problematic because a higher percentage of these popu- 
lations live in areas hosting TRI facilities and they lack 
the salutogenic infrastructure required to act as a buffer 
against exposure to toxic emissions and other environ- 
mental stressors. 

This study had several strengths with the first being that 
all data were from the same year (2010). In studies that 
use census data, this can only be said once a decade. In 
addition, this work can contribute to statewide planning 
as part of PlanMD [68] which is Maryland's comprehen- 
sive plan for sustainable growth and development. While 
these analyses were specific to MD, the methods can be 
used as a template for other states trying to illustrate the 
same relationships among SOD composition, presence of 
pollution-emitting facilities, and health care infrastructure. 
Another notable strength of this study is that it provides 
insight into SOD measures other than race/ethnicity that 
may be a better indicator of spatial disparities in the distri- 
bution of TRI facilities within census tracts. Furthermore, 
our analysis of HPSA versus non-HPSA census tracts 
hosting TRI facilities can be used by the state to guide the 
allocation of resources to help reduce toxic releases in vul- 
nerable communities. 

After completion of this study, there is still a need for 
further analyses. We did not measure actual exposure in 
fenceline communities near each TRI facility or the vari- 
ous buffer zones but rather used distance as a proxy for 
exposure to TRI facilities. In addition, we did not in- 
clude TRI emissions data or measure any health out- 
come data in the buffer zones to show differential 
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health status in conjunction with the location of TRI 
facilities. However, it is worth noting that in terms of ef- 
fect magnitude some findings we identified via multi- 
variate linear regression were not significant, but they 
provide an overall trend in the population under study. 
In reference to HPSA status, additional research is 
needed to address the potential "double disparity" effect 
for non-white and economically disadvantaged commu- 
nities who are living in HPSA tracts that were dispro- 
portionately burdened by TRI facilities and may have 
differential exposure to toxins reported and not re- 
ported by facilities under EPCRA. 

In the future, we plan to incorporate Risk-Screening 
Environmental Indicators (RSEI) and National-Scale Air 
Toxics Assessment (NATA) data from the USEPA in 
order to conduct a more comprehensive analysis. Specif- 
ically, the RSEI database provides information regarding 
the amount of chemical releases, the fate and transport 
of the chemical, the route and extent of human exposure 
to the chemical, the number of people affected, and tox- 
icity [69] which would allow us to better estimate expos- 
ure to TRI emissions. Using NATA data would allow us 
to estimate cancer risk so that we could determine 
whether there is a disparity in cancer risk in communi- 
ties hosting TRI facilities. In addition we plan to further 
examine trends in facility siting and changes in emis- 
sions over time from 1990, 2000, and 2010. 

One final consideration is whether the use of census 
tracts is the most appropriate population measure. While 
there is a wealth of information available at the census 
tract level, they are often not representative of true neigh- 
borhood boundaries. If we decide to represent a real 
neighborhood analysis in the future, we will have to move 
away from the use of zip code and census tract analyses in 
our research particularly for metropolitan statistical areas. 
We see potential in the approach taken by the Baltimore 
Neighborhood Indicators Project [70] to measure and 
track actual' neighborhood level social, environmental, 
and health data. We also believe that the use of planning 
districts or councilmatic districts may be a useful alterna- 
tive to census tracts. 

Conclusion 

This information may be useful to community-based orga- 
nizations seeking to obtain information on the spatial dis- 
tribution of TRI facilities and assistance from federal 
agencies such as the USEPA and the Agency for Toxic 
Substances and Disease Registry (ATSDR) to study the 
negative health impacts of these sites as part of a compre- 
hensive community revitalization program. In addition, 
state agencies such as the MD Department of the Environ- 
ment and the MD DHMH may be able to use the results 
of this study in its efforts to prioritize areas in vulnerable 
communities with a high concentration of TRI facilities 



and toxic releases and leverage state resources to clean up 
areas, improve public health, and enhance quality of life 
and community sustainability. 

Additionally, this work has utility in providing metrics 
for how federal and state regulatory programs are meet- 
ing goals to reduce environmental injustice and environ- 
mental health disparities including cumulative impacts 
of environmental hazards in environmental justice com- 
munities. In addition, this work can contribute to state- 
wide efforts to reduce health disparities and achieve 
health equity through implementation of the Affordable 
Care Act or state laws such as the MD Health Improve- 
ment and Health Disparities Reduction Act. With the 
right investment of ACA resources, hospitalizations, 
emergency room visits, and overall burden of disease re- 
lated to exposure to toxins and other agents could be re- 
duced in overburdened and underserved areas. 

Additional files 



Additional file 1: Table SI. BIC Stepwise Model Selection on Both Directions. 

Additional file 2: Table S2. Confounder Selection (Each possible 
confounder was added sequentially). 
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